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Abstract 

Gradient expansions in quantum transport equations of a Kadanoff-Baym form have been 
reexamined. We have realized that in a consistent approach the expansion should be per- 
formed also inside of the self-energy in the scattering integrals of these equations. In the 
first perturbation order this internal expansion gives new correction terms to the generalized 
Boltzman equation. These correction terms are found here for several typical systems. Possible 
corrections to the theory of a linear response to weak electric fields are also discussed. 
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1 Introduction 



Time dependent transport phenomena in quantum many-body systems can be de- 
scribed by the nonequilibrium Green's function formalism (NGF) of Kadanoff and 
Baym [Q] or Keldysh 0. The Kadanoff-Baym transport equations for nonequilib- 
rium correlation functions can be obtained by analytic continuation to real times 
|| of the Dyson equation for Matsubara Green's functions in purely complex times 
@, H|. The differential form of these equations has been applied in many systems 
B H? H> SI- Usually it is necessary to approximate the equations on several levels 
although in some systems the equations can be directly solved by powerful nu- 
merics ||10||. The integral form of the Kadanoff-Baym equations was less exploited, 



because approximations in a time domain are not so familiar here ||1 1|| . Both these 
approaches have been also used to develop a linear response theory for quantum 
systems in weak dc and ac-electric fields. Here the integral version [12[] seems to be 
more direct than the older differential versions [13|, [14| . 

All functions in the Kadanoff-Baym equations depend separately on two time 
and two space arguments (ri,^;^,^)- This two-argument structure, which result 
in nonequilibrium many-body systems with time or space nonlocal scattering, is the 
main obstacle in solving the transport equations. An approximate one-argument 
form of the equations can be obtained in systems with not very strong interactions 
[T5| . Most older transport methods have a one-argument structure, because they 



usually implicitly consider the presence of weak interactions. 

A different simplification can be obtained in all types of systems if exter- 
nal excitation fields vary slowly in time and space. Then it is useful to subtract 
two equivalent sets of the Kadanoff-Baym equations, differentiated over the first 
(ri,ti) and second variables (t^,^), and perform the so called gradient expan- 
sion m in the new equations transformed to the center-of-mass system (CMS) 
f = (x, X) = (r, t; R, T) = (r x -r 2 ,t 1 -t 2 ; a±a } *i±k) . T n t h e transformed equa- 
tions various terms have the second arguments with big variables (R, T) shifted by 
different fractions of the small variables (r ,t'). A Taylor expansion of these terms 
in the small variables arround the common values of the big variables is the gradient 
expansion. This expansion in powers of derivatives over (R,T) multiplied by (r ,t') 
can be unambiguously and consistently stopped at chosen perturbation orders. In 
this way the two sets of CMS variables (R,T) and (r,t) can be step by step de- 
coupled. By this decoupling memory effects are cuted, so that the nonequilibrium 
dynamics becomes quasi-local and the equations get a quasi-equilibrium form. Far 
from equilibrium this approach evidently fails, because the nonlocal scattering, lead- 
ing (in equilibrium) to quasiparticles with a (k, cj)-dependent self-energy, should be 
reflected in the nonlocality of the nonequilibrium dynamics. 
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In the zeroth order of the gradient expansion the transport equations get a 
local dynamics. Stopping the gradient expansion in the first order gives corrections 
to the equations which partially restore memory effects peculiar to the nonlocal dy- 
namics. These equations are called a generalized Boltzmann equation [Q] (GBE). 
The GBE is not limited to weak interactions, but it can only describe slow dynamics 
close to equilibrium. The weaker is the scattering the better is the description far 
from equilibrium. When the scattering processes can be considered extremally weak, 
then the nonequilibrium correlation functions in the GBE can be substituted by a 
delta-like spectral function multiplied with a distribution function for momenta 
As a result the GBE reduces to the Boltzmann equation (BE), which differs from its 
classical counterpart only by the degeneracy of the described gases. It is interesting 
that the integral version of the quantum transport equations [|l(J approximated by 
low orders of the gradient expansion do not give the BE. It would be also good to 
mention that the gradient expansions in quantum transport equations are analo- 
gous to expansions in the classical Enskog's equation or the more general BBGKY 
equations [|17 |. 

The gradient corrections in the GBE described in the past [|]] do not fully 
reflect the character of scattering processes, because the self-energy in the scattering 
integrals of the transport equations is considered as a structure-less entity. We 
have realized [18|] that consistently performed gradient expansion should include 
also expansion of the self-energy itself, as soon as the self-energy includes additional 
scattering events separated by internal vertices ||]. In the first order this internal 
expansion gives new correction terms to the GBE, which are determined by the 
character of many-body scattering processes. 

In this work phenomenological rules are presented, which allow to perform 
these internal gradient expansions. The rules are applied to consistently derive the 
GBE with all correction terms. These correction terms are found in three examples 
of a self-energy: the averaged T-matrix approximation for a self-energy in electron 
scattering on local potentials, the local T-matrix approximation and the shielded 
potential approximation for a self-energy of interacting spinless Fermions. The 
importance of the internal gradient corrections for a linear response to weak electric 
fields is also briefly discussed. 
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2 Gradient expansions in quantum transport equations 



In slowly changing fields the two differential forms of the Kadanoff-Baym equations 
in QB.7|) can be subtracted and the gradient expansion can be performed in the 
resulting equations. The right sides of the subtracted equations include scattering 
terms of the form 

S a (xi,x 3 ) G /3 (x 3 ,x 2 ) , G a {x u x 3 ) S /3 (x 3 ,x 2 ) , Xi = (r^U) , (1) 

where the analytical structure of the self-energy £ and the Green's function G in 
(P is determined by the index a, /3 = r, a, <, > (see the Appendix A). 

2.1 External expansions 

Consider as an example the gradient expansion in the term Y, a (xi,x 3 ) G^(x 3 ,X2) 
from ([[J), where the many-body structure of S Q is neglected The expansion 
results in the CMS coordinates x — ( x i X) = (r, t; R, T) as follows 

j dx 3 x 3 ) GP(x 3 , x 2 ) = J dx 3 S a (x t - x 3 , Xl — ^ G p f^x 3 - x 2 , X3 ^ X2 ^j 

= J dxY, a (x, X + G? (x - x, X - |) 

= Jdx X) G^x -x,X) + Jdx X) (^y^) ° P ( X ~ X ) 

(-§) + .-, (2) 

The new coordinates are equal to x = X\ — x 2 , x = X\ — x 3 , X = Xl + X2 . A Fourier 
transform of the last expression in (|2|) from the coordinates x = (r,t) to q = (k,u) 
gives the following expansion in a series of Poisson brackets |J 

exp (I D(e,0) G^(0 = G^(0 

z / d£ a (£)dG^(£) 9E°(0^(Q dS a (£)<9G^(Q dE Q (£) ^(Q \ 
+ 2 V dR dk dk dR df du + du dT J + "' 

= £«(£) G*%) + £ G^(0] + - , t=(Q,X) = (k, u; R, T) . (3) 

The gradient expansion of the scattering term H a (xi,x 3 ) G^(x 3 ,x 2 ) has been ob- 
tained in (H) without taking into account the many-body structure of the self-energy 
£ a . Therefore this expansion can be called external. The first order of this expan- 
sion includes only the first Poisson bracket in @. Similarly can be performed the 
gradient expansion in driving terms of the quantum transport equations. 
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2.2 Internal expansions 

In many systems the self-energy is approximated by Feynman diagrams formed by 
ladders or bubbles of Fermion and Boson Green's functions [[l], [|, |19| . This nontrivial 
functional of Green's functions, connected by internal integrations, can substitute 
the self-energy in any step of calculations. The internal structure of the functional 
was not taken into account in the gradient expansion (|3|). Does it contribute by new 
terms in the gradient expansion? It is obvious that this question can be answered, 
if a systematic gradient expansion is performed in the expressions ([I]), where the 
self-energy is substituted by the functional of the Green's functions. 

Consider for simplicity that the self-energy describes electron scattering on the 
localized potentials V(r) = J2n <K r ~ r i)- Assume further that the value Vo is of 
a moderate strength, but the number the random coordinates r« is relatively small. 
Then the suitable self-energy for this problem results by the averaged T-matrix 
approximation |2(J (ATA) 

h) = c{V + V 2 G(tx,t 2 ) + V 3 G(h, t 3 ) G(t 3 , t 2 ) + ...} = c 9(t u t 2 ) . (4) 

Here c represents the weak concentration of the potentials Vo and 9 is the local 
T-matrix. The Green's functions in (|J) depends only on the time variables, because 
the space variables have been integrated out due to the local scattering. 

Analytical continuation of the self-energy (|4]) gives the propagator and corre- 
lation functions for the self-energy in the form (see the Appendix B) 

£ r (ti, h) = c {Vo + V 2 G r (t u t 2 ) + V 3 G r (h, t 3 ) G r (t 3 , t 2 ) + ...} = c 9 r {t u t 2 ) , (5) 
£<(*!, t 2 ) = c Pfah) G<(t 3 , U) 9 a (U, t 2 ) = c 9<{t 1 ,t 2 ) . (6) 

Each term from the expressions for E r,< in (^-^j is formed by several propagators 
or a correlation function with time arguments in a 'series' (see the Appendix B). 

Let us study in details the gradient expansion of the scattering terms in ([I]) 
for the self-energy in (0-|6|). Assume for the beginning that the self-energy in the 
first expression from (|TJ) is substituted by one of the terms from (|5|-§D, which has 
two Green's functions A£(£i,t 3 ) = c Vq 3 G a (tx,i 2 ) G^fo,^) (no index is used at 
the contribution to the self-energy AS to show its analytical structure). Then 
the expression AS(t!,f 2 ) C 7 (t 2 ,^3) can be transformed to the CMS coordinates as 
follows (we suppress the prefactor c V 3 and neglect the fact that the function G 7 
depends also on space variables) 

G a (h,i 2 ) GP(i 2 ,t 3 ) G\i 3 ,U) 
= G a (t, - i 2 , *±*l) G 13 (t 2 - i 3 , ^) G J (t 3 - U , ^) 
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T-Ti-Tu , T/A „ fl /_ _ , T-Tj-Tu Tf 



G a ( 77 , T + — — + -jf- ) G? f H ,T + 



2 J V 2 2 / 

x C 7 - f/ - f n , T - y - yj , a, /?, 7 = r, a, <, > , (7) 

where it holds (r, T) = (£i — £ 4 , ^j ± ) and fi = t\ — i 2 , f n = i 2 — £3. In the 
expression (0) a Taylor expansion in the r-coordinates can be performed around the 
T-coordinates as in ([|). In a first order each of the Green's functions differentiated 
over T has multiplicative coefficients formed by the r-coordinates of the remaining 
two Green's functions in a series. Therefore in this structure of arguments the 
Green's functions are equivalent from the point of view of gradient expansions. 

After a Fourier transform over the small variables a symmetrical expression 
in the three functions G a,l3r( with respect to the derivatives ^ can be obtained. 
Therefore the complete gradient expansion up to the linear order can be shortly 
written as follows (the space variables have been included in G 1 ) 

cV 3 [G a (u,T) G%co,T) G\i)+ l - [G a ^,T) G?(u,T) ,G 7 (0] 

+ l - T), T)} G 7 (0} , £ = (k, u- R, T) , (8) 

where the Poisson brackets have been used. Since the space coordinates (k, R) are 
integrated out in the present self-energy, gradient expansions cannot be performed 
in these variables. But in a general case both pairs of coordinates contribute in the 
way shown in (|8|). 

The second term in (§) resulted as in (|3|) by the (external) gradient expan- 
sion of the expression A£(£i,t 2 ) G 7 ^,^)- The last term in (|j) resulted by the 
gradient expansion of the internal structure of the self-energy contribution AS. 
Therefore this expansion and the resulting term can be called internal. In (§) only 
contributions to the self-energy diagrams with two Green's functions were included. 
Analogously can be performed the gradient expansion in terms with any number of 
Green's functions. This problem is solved in the next section in details. 

The internal gradient corrections can be easily physically understood on the 
previous example, where scattering processes can be seen as many consequent events 
on the same center. If the system is excited by a time dependent field, then scattering 
conditions on the center can change between these consequent scattering events. 
Corrections to these changed scattering conditions are represented by the internal 
correction terms. For centers which are little smeared in space, internal gradient 
corrections result nonzero even in excitation by static fields. 
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3 Gradient expansions in general 



A gradient expansion in quantum transport equations can be unambiguously repre- 
sented by a series of Poisson brackets of an increasing order n. Physically this is a 
consistent expansion in powers of space (time) inhomogeneities, since the n-th order 
Poisson brackets could be appreciated by terms of the form (ko~) n , where k is the 
inverse mean free path and a is the range of space inhomogeneity (analogous terms 
apply for the time inhomogeneity). If ka <C 1 then the expansion can be stopped in 
the first order (n = 1), which is equivalent to inclusion of both external and internal 
corrections in the GBE (see also (||), @ and fl32|)). 

Gradient expansions can be performed also in classical transport equations 
describing dense systems. Such systems have been firstly approximately studied by 
the Enskog's equation ||T7|| , which generalizes the Boltzman equation by taking finite 
volumes of scattering particles. Since only binary collisions are included here, like in 
the BE, gradient expansions in this equation are analogous to the external expansions 
in quantum transport equations, giving the uncomplete quantum GBE [|TJ. Later 
on classical dense systems have been described by the so called BBGKY hierarchy 



of kinetic equations 17], which can include multiple encounters of particles. Under 



some assumptions this set of equations can be reduced to a classical generalized 



Boltzman equation |fL7|| , which has scattering integrals with a structure analogous 
to that in the quantum transport equations. Therefore gradient expansions in this 
classical GBE are analogous to both the external and internal expansions leading 
to the complete quantum GBE (in the classical GBE gradient expansions can be 
performed, while the quantum GBE is the result of the gradient expansions). 

3.1 General rules for gradient expansions 

We formulate a set of phenomenological rules for performing external and internal 
gradient expansions in terms like in (|I]), appearing in quantum transport equations. 
From now the gradient expansions are stopped after the first perturbation order 
(higher order contributions can be found analogously). The following qualitatively 
different points specify and summarize the necessary steps for gradient expansions 
in the concret terms: 

1. Analytic continuation to real times of the term is performed, to get an expres- 
sion formed by propagators and correlation functions. 

2. Each of the self-energy functions Y7' a ' < > > is resolved into a functional of full 
Fermion and Boson Green's functions G r,a,<,> and undressed matrix elements. 

3. The whole term is transformed into CMS coordinates, and the big coordinates 
X are linearized in the conjugated small coordinates x. After a Fourier trans- 
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form over the coordinates x, the linearization prefactors x become derivatives 
d q , which produce a series of new terms. In each of these new terms just two 
derivatives appear (over X and q). The expansion should be performed on a 
hierarchical structure of levels, going more and more inside to the structure of 
E (the vertices are the landmarks). Corrections for higher levels are done on 
lower levels. All objects, which depend on some of the coordinates, must be 
differentiated. The lowest objects are the full Green's functions and undressed 
coupling matrix elements. 

4. 'Parallel' objects with the same coordinates (x;, Xj), but the order of 
considered in the derivatives as a single differentiated object. At the lower level 
some of the objects with the coordinates (xi,Xj) can have still other internal 
coordinates, which can give further gradient corrections. 

5. 'Serial' objects with arguments like in (0) are differentiated in such a way that 
differentiation of one object over a big coordinate X is accompanied by differ- 
entiation of the other objects in a series (one by one) over the conjugated small 
coordinates q. Sign prefactors depend on the sequence of objects. Realization 
of this rule for terms with many Green's functions in a series can be done by 
a second functional derivative of the term over the objects in the series. This 
derivative is multiplied by a Poisson bracket of the two differentiating objects 
in the functional derivative. 

In the previous section these rules have been already implicitly applied on a 
simple algebraic term from a self-energy (see (§)). Analogously can be dealed other 
such algebraic terms or complex recursive terms. Before we come to these expansions 
let us still present one theorem. Application of the above rules to scattering terms 
from transport equations, which include a self-energy as in gives the following 
formal expansion (stopped in the first perturbation order) 

G a (0 S(0 + \ [G°(0, + G a (0 F,[E](0 . (9) 

Here [£](£) represents the internal gradient expansion in the self-energy E. Since 
the form @ is fully general, it can be taken as a theorem: 

• Internal expansion of a term with a self-energy results by the substitution of 
^(0 by £(£) + i*i [£](£) in the zeroth order term. 

From this theorem it follows that the internal gradient expansion in terms with a 
self-energy can be directly found from the function Fj[E](£). 

3.2 Expansions of a self— energy 

We can concentrate on this function Fj [£](£) and evaluate it for several typical 
examples of a self-energy. 
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3.2.1 Static averaged T matrix approximation 

The first example concerns the ATA self-energy (f|) for electron scattering on local 
potentials. After application of the rules 1-3. the zeroth order propagator and 
correlation parts for the self-energy in a frequency representation results 

*w= c T=<Ffawr enu ' T >> (10) 

£<(cu,T) = c9 r {u,T) G<(cj,T) 6 a {u,T) = c 9 < (u,T) . (11) 
Application of the rule 5. to the retarded part in (|5]) gives the first order term 

m r m = ^£f(0 \ [^(0,^(01 = o , (12) 

since only functional derivatives of Y7 over the propagator Green's functions G r can 
be applied here. Because the Poisson bracket from equivalent objects is zero, the 
function Fj[£ r ](£) does not contribute to the internal expansion. Analogously can 
be found the internal term for the correlated part in ([]) 

+ ^^ t 2^^ = - = 
= -clm((6 r (uj,T)) 2 e a (uj,T) [GT(u, T), G<(u, T)}) - °- \9 r (co, T)| 4 G<(u,T) 
r duo 1 ( dA(uj,T)dA(u,T) dA(u,T) dA(u,T) \ 

The term i ? i [E > ](a', T) can be evaluated in the same way. In the first expression from 
(P^) only the nonzero functional derivatives have been considered, which result by 
differentiation over Green's functions of different analytical structures. The order of 
derivatives and terms in Poisson brackets is the same as the order of these functions 
in (|B|). The second expression results by per-partes integration in the Poisson bracket 
[G r (£)> G a (£)] and some simple algebra. 

3.2.2 Dynamic T matrix approximation 

Another T-matrix approximation is used ||1], if mutual interaction of electrons 
are studied. This approximation gives a self-energy of a very similar structure to 
the ATA self-energy in (fD, but internal dynamics is more complicated here. For 
spinless Fermions the self-energy can be written as folows |J 

E(ti,t 2 ) = -«e(t 1 ,t 2 )G(t 2 ,t 1 ) , (14) 
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where is the dynamic T-matrix 

e(ti, t 2 ) = v 5(h -t 2 ) + v R (t u i 3 ) e(t 3 , h) , 

R {hM)=iG{t x ,t 2 ) G(t u t 2 ) . (15) 

In this self-energy the singular Hartree term is included (t 2 — > t x ), but the exchange 
terms have been neglected for simplicity. 

Analytical continuation to real times of QT5| ) gives the propagators and corre- 
lation functions (see the Appendix B) 

Ro(ti, h) = G r (ti, t 2 ) G (t%, t 2 ) — G < (ti, t 2 ) G r (ti,t 2 ) , 

R^(t 1 ,t 2 ) = G < (t 1 ,t 2 )G < (t 1 ,t 2 ) (16) 



and 



Q r {h, t 2 ) = v s(h -t 2 ) + v R r (t u i 3 ) e r (t 3 , t 2 ) , 

Q < (t 1 ,t 2 ) = V[Q r (t 1 ,i 3 ) R<(i 3 ,t 2 ) + Q < (t 1 ,i 3 ) R a (i 3 ,t 2 )} • (17) 



After application of the rule 3. in (16-17) the zeroth order terms can be completed 



in a form similar to the ^-functions in (|T0|-|TT|) 



e>,T)= J , 0<(^T) = e>,T) R<(co,T) Q a (u,T) . (18) 

1 - Rq[uj, 1 ) V 

The zeroth order contributions to the propagator and correlated parts of the self- 
energy (|T1| ) results from these functions as follows 

£ r (cu, T) = &'{uj + to, T) G < (l>, T) - Q < {oj + uj, T) G a {u, T) , 

S<(o;,T) = Q<(uj + 0,T) G > (u,T) . (19) 

We can continue with the internal gradient expansion in the self-energy flU]). 
The separate Green's function in (19) is on the highest level, so it is excluded from 



the internal expansion. The remaining functions r and < could be independently 
expanded after the rule 4. and give Fj[O r ](£) and F i [6 < ](^). 



Time arguments of the objects Ro(ti,t 2 ) in the function Q(ti,t 2 ) from (jTJD 
are in a series, so that the rule 5. can be directly applied. The expansion of the 
propagator function in terms of R r results zero, similarly as in (|T^) 

Fi[er]{0 = ^w {0 \ = o • (20) 
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The internal expansions of the correlation function B< is formed by terms analogous 
those in fO) 



- Im((6 r (a;,T)) 2 e a (a;,r) [^(w, T), J?<( W , T)]) - 2 |6>, T)| 4 i£(w,T) 



rift 1 / dlmflg (g » r ) 9Im^ (a; , T) glmgS (g , T) dIm#S (a; , T) ' 



(21) 



The function i ? i [B > (^)] results by the change of the index < by > in all places of 
Finally the internal expansions -Fj[£ r ](£) and F i [S < ](^) can be obtained, if the 
functions r (£) and < (O in the zeroth order self-energy ( Jl9| ) are substituted by 
F,[9 r ](e) and F t [Q<](0 from (^-0) 



F 4 [E r ] (c, T) = F t [Q<] (u + u, T) G a (a>, T) , 

i^E^T) = F l [e<}(u + cu,T) G>(u,T) . (22) 

Here the fact that Fj[O r ](£) = in fl2H|) has been taken into account. 
3.2.3 Shielded potential approximation 

Particle interactions are often long range, like in the Coulomb potential. If free 
particles are available in the system, then screening can shorten the range of the 
interaction (give a much more localized potential). The shielded potential approxi- 
mation of a self-energy ^ does not include ladders, like in the above studied exam- 
ples, but rows of electron-hole bubbles screening the potential. Localization of the 
potential by screening has a dynamical character. Therefore it would be demanding 
to know time dependent behavior of screening in nonequilibrium processes. Re- 



cently dynamics of ultrafast screening processes have been studied by NGF JH], |22 
Close to equilibrium these systems can be described by the GBE, where the internal 
corrections to the screening dynamics can play an important role. 

The shielded approximation for the self-energy can be written in the form |I[ 

S(l,2)=zK(l,2)G(l,2), (23) 

where the singular Fock term is included (t% — > tf). In the lowest order of a pertur- 
bation theory, the screened potential V s is related to the unscreened one V by the 
polarization function L Q as follows 
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V S {1, 2) = V(l, 2) - V (1, 3) L (3, 4) V s (l, 2) , V(l, 2) = F(n - r 2 ) 5{U - t 2 ) , 

Lo(l,2)=iG(l,2)G(2,l) ■ (24) 

The propagators and correlation functions for the polarization function L in 
(|24l) can be found by the rules in the Appendix B 

U (l, 2) = -((7(1, 2) G<(2, 1) + G<(2, 1) G a (2, 1)) , 

L<(1,2) = -G < (1,2)G > (2,1), (25) 

and the functions VJ and result analogously as r and < in (|i~7l). 

Therefore after application of the rule 3. in Q2"5D, the zeroth order of the gradi- 
ent expansion can be obtained as in (p~8|) (the full CMS coordinates £ = (k,u; R,T) 
are considered, since the interaction is nonlocal) 

Ki0 = i + m)v(k) ' Vs<{ ° = v:{0 Lo<(0 v:{0 ■ (26) 

Similarly the zeroth order gradient contributions to the propagator and correlated 
parts of the self-energy (p3|) result 



Z r (k,u;R,T) = G r (k-k,u-u;R,T) V?(k, w; R,T) 
-G < (k-k,u-0;R,T) V s a (k,u; R,T) , 

S< (fc, u\ R, T) = G< (Jfc - k, uj - oo; R, T) V< (k, u- R, T) . (27) 



The structure of the internal corrections to the self-energy (|24|) is also similar 
to the previous example. The main difference is in that here the potential V(k) 
depends on the wave vector k, so that new corrections terms are included 

+ Mv (0 \ imhV{k)] + wh {& 5 m)] = ° ' (28) 

The potential V in the last two terms is either more 'right' or more 'left' than L in 



the series ( p4[) (see also comment below (p~3|) ) . Since the second functional derivative 



is symmetrical in the differentiating functions, these two terms diminish each other. 
In the Poisson brackets from fl2"8p derivatives over space-momentum coordinates are 
also performed (see (§)), but the potential V(k) can be differentiated only over the 
momentum. 
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The correlation function V< gives more interesting internal contributions. We 
write only terms with nonzero Poisson brackets and also neglect terms analogous 
those in (|28f) , where both functions originate from one of the propagators Vp a . The 
fact that the differentiating functions V, Lq 0, originate from different full potential 
propagators Vp a is symbolized by a small index at the potentials V^ r \ V^ a \ which 
show the side they originate from. The nonzero terms are 

Fi[Vs<m = ^< (0 ^ [zmL ° <(0] + t [L ° {0 > Lm 

S 2 V< _ i 



5Ll5L a Q J 2 
S 2 V< _ % 5 2 V< 



+ 6$fa>® \ ^ V ^ + 5^k^ \ ^ L ^\ • (29) 



The first three terms in the right side of fl2"P| ) are analogous those in fl2"Ip, but the 
additional space-momentum derivatives present in the Poisson brackets here. The 
next two terms in fl2D|) are evidently nonzero, because the derivatives and the Poisson 
brackets are different and nonzero. In the remaining two terms the Poisson brackets 
are opposite each other, but the functional derivatives are different, because the 
potentials are taken from different sides. 

We can evaluate the above functional derivatives and collect all the terms 

mm) = i^{(K r (e)) 2 K a (C) {[L r (O,m)} - y^yi [v(k),L<(oM 

(30) 

The function -^[^(01 results by the change of the index < by > in all places of 
©■ 

As before the internal expansions Fj[£ r ](£) and Fj[E < ](^) can be obtained, if 
the functions VJ' a (£) and ^ <5> (0 in the zeroth order self-energy ( p7|) are substituted 
by Fi[V s r ' a ](0 and &[¥<>>](£) from (HfOD 



Fi[J: r ](k,uj;R,T) = G r (k-k,u-u;R,T) Fi[V>](k, cu; R, T) , 
F i [S<](A:,aj;i?,T) = G < {k -k,u - oo; R,T) Fi[V<](k,u; R,T) . (31) 
The fact that -Fi[V7' a ](£) = in (p8|) has been again taken into account. 



13 



4 Transport equations with internal corrections 



The internal corrections might appear in all kind of equations derived by the gradient 
expansions from the Kadanoff-Baym equations. 

4.1 Generalized Boltzmann equation 

The complete generalized Boltzman equation with both the external and internal 
gradient corrections can be found directly, when the above theorem is implemented 
in its derivation from the Appendix B. The theorem says that the internal terms 
result by the substitution of £ <(>) (0 by £ <(>) (£) +F i [£< > )](£) in the zeroth order 
terms. Therefore after a Fourier transform the complete GBE gets the form 

( 9 . P'Vfl jr (ft vv | 9U eff (Q 0\ 

-[ReS r (0,G<(0] + [ReG"-(0,S<(0] 

= - (£>(£) + F i [E>(0]) G<(0 + (E<(fl + fi[E<(fl]) G>(0 • (32) 

The equation for the correlation function G > (£) can be found similarly. The singular 
(Hartree-Fock) contributions in ([32]) are included in the effective potential C4//(0- 
The propag ator functions G r ^(£) and £ r(a) (£) 

in (|3*^) are related to the correlated 
parts G<(£), G>(£) and £<(£), S > (0 by the Hilbert transform (Q. Therefore it 



is not necessary to find separate transport equations for these propagator functions, 
which would eventually include the nonzero internal terms Fj[E r ( a )(£)]. The internal 
terms in fl32|) contribute to the dynamics (not the renormalization) of the studied 
system, so that they cannot be neglected in the transport equations. Conservation 
laws for the GBE in ([32]) could be proven similarly as for the exact Kadanoff-Baym 
equations [|TJ. 

In the above example of the ATA self-energy, the expressions Fj[£ <(> )(£;)] in 
(p^) should be substituted by (|T3|). For the dynamic T-matrix approximation the 
terms (p2]) fulfill this role (the Hartree term should be taken once). In both these 
examples the CMS coordinates in the functions Fj[S < ^ > )(c<j, T)\ are reduced with 



respect to the other terms in (|32]) , where it is £ = (k, oj; R, T). This reduction results 
from the local form of the T-matrix approximations studied here. In the shielded 
potential approximation for the self-energy, the expressions -Fj[S <( ^ > )(^)] from (|3"lD 
with full CMS coordinates can be used in the GBE. In all these cases the complete 
GBE with internal correction terms results quite complicated. Nevertheless, we 
believe that it can be handled by some approximate numerical methods. 



14 



4.2 Linearized transport equations 

The internal gradient corrections can appear also in linearized transport equations 



derived from the GBE in weak cfc-electric fields [fL3 |. It has been already men- 
tioned p3| that these linearized equations should include further correction terms 
for complicated scattering. 

The left side of the GBE in (^) looks in dc-electric fields as follows 

(^ + > +erE »>'( v * + eE °S) G<(0 ' (33) 

To diminish the explicit time dependence in ([$3|), not present in dissipative systems 
in weak cfc-fields, the following transform should be performed |L3[ 



Q^k + eEoT, _^_ + eE -V Q , (34) 

where E is the intensity of the dc-electric field. Application of this transform to 
the Poisson brackets in the second line of (^) gives mixed terms in the following 

way 

[A, B] -> [A, B] + eE • V Q B - V Q A ^) . (35) 

If this transform is applied also to the Poisson brackets in the internal corrections 
F i [E < (»(0] in (|32|), further new terms can result. 

When the self-energy E(£) depends only on (a;, T) variables, like in the space 
localized scattering ([|) or (|T4|), then only the external terms dRe g w ^ 9G d ^ k ^ and 
9T:< du}^ dReGJ^) ; result by the transform (^). Application of this transform in 
the terms Fj from (|T3"D, (^) give nothing, because only fc-independent 

Green's functions are present there. If a nonlocal ATA self-energy is used, 
generalizing (for neutral smeared imperfections) the local form ([|), then new terms 
would result from application of the transform (|35| ) in the A;-dependent internal 
corrections Fi[E < ^ > - ) (^)]. The nonlocal scattering results for example also on charged 
impurities , where it is reasonable to screen the impurity potential , similarly 



as in the self-energy (|2B|). In both these examples internal and external corrections 
in (|3"2| ) are nonzero. Since evaluation of these terms is direct, we do not write them 
as well as the resulting complicated linearized equations. 

The question is how the internal corrections can contribute in the case of 
the linear response to weak ac-electric fields, where we would expect that also the 
space-local interactions (|4]), ( |T4D give new correction terms. The linearized trans- 



port equations in the ac-electric fields have been studied W%\, but the transform to 
new coordinates ]TB[ was not performed. Therefore no gradient corrections would 
seemingly contribute by new terms. In fact it is probably hard or not unique to 
find the above transform in the ac-case. Moreover in the ac-case it is not sufficient 
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to stop the gradient expansions in the lowest orders. Therefore a new gauge invari- 
ant approach has been developed f27j to study linear response to the weak ac and 



cfc-electric fields, where no additional transforms are necessary. Unfortunately, the 
resulting equations are still quite complicated. Recently a relatively simple consis- 
tent approach has been devised jn|, which starts from the integral version of the 
Kadanoff-Baym equations. 

5 Conclusion 

We have found new gradient corrections in the generalized Boltzmann equation 
These corrections result if the gradient expansion is performed also inside of the 
self-energy in scattering integrals of the quantum transport equations. We call 
these corrections internal, because they reflect the many-body character of scatter- 
ing processes represented by the internal structure of the self-energy. Analogously 
the standard gradient corrections [||], which do not take into account the internal 
structure of the self-energy, are called here external. 

The generalized Boltzmann equation with all correction terms has been de- 
rived. The internal corrections to the GBE have been calculated for electron scat- 
tering on localized static potentials, which is described by the ATA self-energy. More 
complex corrections have been obtained for interacting spinless Fermions, where the 
self-energy is described either by a local T-matrix approximation or by a nonlocal 
shielded potential approximation. We believe that the GBE with the new correction 
terms might be a proper tool for studies of relaxation to equilibrium in systems with 
nontrivial electron interactions. We are planing to investigate some of these systems 
in future. 

We have also discussed the importance of the internal corrections in the lin- 



earized transport equations in weak electric fields [[L3|, which can be derived from 
the GBE. The presence of new correction terms in the cfc-version of these equations 
has been clarified in the above studied examples. The internal corrections might be 
important in many other physical problems, where the self-energy includes multiple 
scattering events. 
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Appendix A 

The causal Fermion (O = ip) or Boson (O = A) Green's functions in real times are 
defined by j|, [5] (Matsubara Green's functions in complex times look analogously) 

G\l,2) = -i < T[0(1) Ot(2)] > , j = (r v t,) , (j = 1, 2) . (A.l) 
Correlation functions are related to the causal function as follows 
ih G*(l,2) = G > (1,2) =< 0(1) O f (2) > , h > t 2 , 
=F i% 2) = G < (1, 2) =< O f (2) 0(1) > , *i < t 2 , (A.2) 

where the upper (lower) sign applies to Fermions (Bosons). 

The retarded and advanced Green's functions are defined by 

G r (l,2) = -ifl(l-2) [G>(1,2)±G<(1,2)] , 
n 

G\l,2)= l - 6(2 - 1) [G>(1, 2) ± G<(1, 2)] , (A.3) 

where the theta function is 6{t) = 0, t < 0; 6>(t) = 1, t > 0. 

In equilibrium and space homogeneous systems the Green's functions depend 
only on the difference of coordinates (r, t) = (r\ — r 2 ,£i —t 2 ), so that they can be 
easily Fourier transformed to the (/c, a>) -represent at ion as follows 

G(k, w) = f d n r J dt exp(i(ut - r ■ k)) G(n - r 2 ; t x - t 2 ) . (A.4) 

Then the Fermion and Boson correlation functions can be expressed as jl| 

G < (k,u) = n FjB (huj) A(k,u) , G y (k, u) = (1 =f n F>B {huj)) A(k, uj) , (A. 5) 

where n F , n B denote the Fermi-Dirac and Bose-Einstein distributions 

n FyB (hu) = ^} 

ekT ± l 

and the spectral function is defined by 

A(k,u) = -2 Im G r {k,u) = G > {k,u) ± G < {k,u) . (A.6) 

The retarded Green's function can be calculated from the spectral function 
as follows ( ( |A.6| - ATT ) hold also in full CMS coordinates £ = (k, u; R, T)) 



J-oo 2ir u — uj + id 



Similar formulas can be applied also for the self-energy. 
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Appendix B 

The nonequilibrium Green's functions can be found by analytical continuation to 
real times of the Matsubara Green's functions in complex times |1|]. In NGF it 
is often necessary to find the propagator or correlation part of a combination of 
functions. An example is the product 

A(l,2) = 3(1,2)^(1,2), (B.l) 

where A, B, C are one-particle causal Green's functions or self-energies. The 
required functions can be found by LW rules ||, where the signs and prefactors 
result from the definitions ( |A.2| - |A.3| ). We have found the expressions for ( |B.1| ) in 
two cases, where the functions A, B, C correspond either to Fermions (F) or to 
Bosons (B) as follows: (1) (A, B, C) = (F, F, B) (an example is an electron-phonon 
self-energy ||) or (2) (A, B,C) = (B, F, F) (an example is the function R in (|T5|)). 
In both these cases the expressions result in the form 

A<(1,2) = -i B<(1,2) C<(1,2) , A>(1,2) = -i B>(1,2) C>(1,2) , 

A r (l,2) = -i (iT(l,2) C>(1,2) - £<(1,2) C r (l,2)) , 

A a {l, 2) = -% (B a (l, 2) C>(1, 2) - B < (1, 2) C a (l, 2)) . (B.2) 

Similarly can be found the propagators and correlation functions for the expression 

A(1,2) = B(1,2) C(2,l) , (B.3) 

where the following possibilities have been chosen: (1) (A,B,C) = (F,B,F) (an 
example is the self-energy Ql4]) ) or (2) (A,B,C) = (B,F,F) (an example is the 
electron-hole bubble L in (0)). These possibilities, which, but the order of B, F 
in (1), are the same as in ( B.2|) , give the identities ((1) for (— ), (2) for (+)) 

A<(1,2) = i S<(1,2) C>(2,1) , A>(1,2) = i B>(1,2) C<(2,1) , 

A r (l, 2) = i (B r (l, 2) C<(2, 1) T B<(1, 2) C a (2, 1)) , 

A a {l,2) = t {B a {l,2) C < {2,1)tB < {1,2) C r {2,l)) . (B.4) 

The following structure appears also in most formulas 

A(l,2) = £(1,3) C(3,2) , (B.5) 

where the bars over the arguments mean integration over the whole real axis of these 
coordinates. The expressions result the same, irrespective of the types of involved 
functions 

A<(1, 2) = B r (l, 3) C<(3, 2) + B < (1, 3) C a (3, 2) , 
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A>(1, 2) = B r (l, 3) C>(3, 2) + B>(1, 3) C a (3, 2) , 
A r (l, 2) = B r {l, 3) C r (3, 2) , A a (l, 2) = £ a (l, 3) C a (3, 2) . (B.6) 



In the text we use the term 'parallel' for the structures of arguments in ( p.lHB.4|) , 
while the structure ( |B.5HB.6| ) is termed 'serial'. 



The Kadanoff-Baym equations can be found by application of the rules ( |B.6| ) 
to the differential Dyson equation, which can be written in two forms. If we take 
into account that (Gq 1 ) < (1,2) = and 5 < (1 — 2) = 0, then the two forms of the 
Kadanoff-Baym equations look as follows 

(G r )-\l,3) G<(3,2) = 57(1,3) G<(3, 2) + S<(1, 3) G a (3,2) , 

G<(1, 3) (Gg)" 1 ^, 2) = (7(1, 3) S<(3, 2) + G<(1, 3) 57(3, 2) . (B.7) 

Analogous equations can be obtained for the correlation function G > . 

The generalized Boltzmann equation jl] can be derived by subtraction of the 
two sets of equations (|B.7|) . In the resulting quantum transport equations for G K,> 



it is necessary to introduce the CMS coordinates and perform the gradient expan- 
sion up to the first order. Then a Fourier transform over the small coordinates is 
performed. It is helpful to resolve the propagators G r,a , 57' a from the right side of 
( |B.7| ) into the real and imaginary parts. The imaginary parts of these propagators 
can be resolved with the help of the identity ( |A.6| ) (in the CMS coordinates). Then 



the terms with equal correlation signs <, < and >, > fall out from the scattering side 
of the new equations and the GBE easily results (see its completed form in (|3"2])). 
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